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Abstract 

A recent result on size functions is extended to higher homology modules: the persistent 
homology based on a multidimensional measuring function is reduced to a 1-dimensional 
one. This leads to a stable distance for multidimensional persistent homology. Some reflec- 
tions on i-essentiality of homological critical values conclude the paper. 
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1 Introduction 

Topological Persistence started ante litteram in 1991 with P. Frosini, who introduced the concept 
of Size Function [I1],[I71 Sect. 8.4], a topological-geometrical tool for describing, analyzing and 
comparing shapes. This was actually the origin of rather large experimental research (|19[ 
121)). Size functions were generalized by the same School in two directions: Size Homotopy 
Groups (already in a multidimensional setting!) and Size Functor [5]. 
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At about the same time, Persistent Homology was independently introduced [101 E] (see 
also [51[S]). All these theories have substantially the same target: shape recognition. They are 
constructed on some topological features of lower level sets of a continuous real-valued function 
defined on the object of interest. They also share an important advantage with respect to 
other methods of pattern recognition: they capture qualitative aspects of shape in a formal 
quantitative way; so, they turn out to be particularly suited to the analysis of "natural" shapes 
(blood cells, signatures, gestures, melanocytic lesions, ...). Retrospectively, a size function 
is identifiable with the rank of a 0-th persistent homology module, while the first persistent 
homology module is the Abelianization of the first size homotopy group [14J, and the size functor 
[2] is a functorial formalization of the direct sum of persistent homology modules. 

The results obtained recently, involving the construction of size functions related to mul- 
tidimensional measuring functions, lead us to the same generalization to persistent homology 
modules, which is the goal of this paper. As far as Size Theory is concerned, the main reason for 
such a generalization is that there are shape features, that have a multidimensional nature (such 
as color) and whose description can be done necessarily by a multidimensional measuring func- 
tion. Moreover, there are shapes, which cannot be discriminated by n size functions related to n 
different real-valued measuring functions, but can be distinguished by the size function related 
to the n-dimensional measuring function of which those are the components (see Section [5|). As 
mentioned in [151 Section 2.5], the study of multidimensional persistence has strong motivations, 
but some objective obstacles. This paper wants to pave a way out of these difficulties. 

After recalling some basic notions about multidimensional size functions and 1-dimensional 
persistent homology in Section [21 we adapt the arguments of 1 to multidimensional persistent 
homology in Section [3l for proving our main result (Theorem [2]) . This is a reduction theorem, 
which takes the detection of discontinuity points back to the case of 1-dimensional persistent 
homology. This seems to overcome the pessimistic final considerations of [3, Section 6] on the 
structure of the functions px,i- In fact, although the sets, on which the functions are constant, 
are much more complicated than the triangles typical of the 1-dimensional case, they reduce 
to them when properly "sliced" by a suitable foliation. Stable distances on the leaves of the 
foliation define (and approximate) a global distance for rank invariants. Examples and further 
remarks on a different kind of reduction conclude the paper. 
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2 Basic notions 



In the first part of this section we'll recall briefly the concept of multidimensional size functions 
and we'll state the theorem that gives us the tools to calculate them (Theorem [1]). It asserts, 
indeed, that a suitable planes' foliation of a 2n-dimensional real space makes an n-dimensional 
size function equal to a 1-dimensional in correspondence of each plane [1]. In the second part 
we shall review the definitions of persistent homology module and related concepts [5j . 

2.1 Multidimensional Size P\inctions and 1-dimensional reduction 

In Multidimensional Size Theory, any pair {X, /) , where X is a non-empty compact and locally 
connected Hausdorff space, and / = (/i, . . . , /„) : X K" is a continuous function, is called a 
size pair. The function / is called an n-dimensional measuring function. The following relations 
^ and -< are defined in R": for u = (ui, . . . , u„) and v ~ (ui, . . . , w„), we say u ^ v (resp. 
u < v) ii and only if Uj < Vj (resp. Uj < Vj) for every index j = 1, . . . ,n. For every n-tuple 
u = (ui, . . . , u„) e R", let X{f ^ u) be the set {P e X : fj{P) < , j = 1, . . . , n} and let A+ 
be the open set {{u, u) e R" x R" : u w}. 

Definition 1. For every n-tuple v — (wi, . . . , w„) G R", we say that two points P,QCzX are 
if ^ t?) -connected if and only if a connected subset of X{f ^ v) exists, containing P and Q. 

Definition 2. The (n-dimensional) size function associated with the size pair (X, f ) is the 
function /) • ^ -f*^' defined by setting £^-^ j<^(u,v) equal to the number of equivalence 
classes in which the set X{f ^ u) is divided by the (/ ^ if) -connectedness relation. 

An analogous definition for multidimensional persistent homology will be given in Definition 

m 

The main goal of T for size functions, and of the present paper for persistent homology, is 
to reduce computation from the multidimensional to the 1-dimensional case. This is possible 
through particular foliations of R" by half-planes. They are determined by what are called 
"admissible" vector pairs. 

Definition 3. For every unit vector I — (Ji, . . . , 1^) in R" such that Ij > for j — 1, . . . , n, 

and for every vector b — (6i, . . . , 5„) in R" such that ^ bj — 0, we shall say that the pair {I, b) 
is admissible. We shall denote the set of all admissible pairs in R" x R" by Adrrin. Given 
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an admissible pair {l,b), we define the half-plane ■""(fg) 'in K" x R" by the following parametric 
equations: 

u = si + b 
V = tl + b 
for s, i S R, with s < t. 

The motivation for the previous definition is the fact that for every (u, u) G A+ there exists 
exactly one admissible pair (l,b) such that (w, w) G ''''(fg) )Xi Pi'op.l]. The following Lemma is 
substantially contained in the proof of 1, Thm. 3]. 

Lemma 1. Let (I, b) be an admissible pair and g : X ^ W be defined by setting 

'fAP)-bj 



g{P) ^ max 

j = l,...,n Ij 

Then, for every (u, v) = (si + b, tl + b) G '^(^ihy following equalities hold: 

X{f^u)^X{g<s), X{f<v)=X{g<t) 
Proof. For every u = (ui, . . . , u„) G R", with Uj = slj + bj,j — 1, . . . , n, it holds that 

X{f^u) = {P^X:f,{P)<Uj,j = l,...,n} 

= {PeX ■.fjiP)<sl,+b„j^l,...,n} 

= {PGX:M£l_^<,,,^i,...,„} 

h 

= X{g<s) 

Analogously, for every H = [vi,.. . G R", with Vj — tlj +bj,j = 1, . . . ,n, it holds that 
X{f^^^X{g<t). □ 

From that, there follows the main theorem of [1]: 

Theorem 1. Let {l,b) and g be defined as in Lemma{^ Then the equality 

holds for every (u, H) — (si + b,tl + b) G '''(f • 
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This is indeed the theorem that we are going to extend, in Section[31 to persistent homology of 
aU degrees. Its importance resides in the fact that essential discontinuity points ( "cornerpoints" 
in the terminology of Size Theory) are the key to a stable distance between size functions. Unfor- 
tunately, cornerpoints do not form, in general, discrete sets in the multidimensional case. This 
theorem makes it possible to find them "slice by slice" with the familiar technique of dimension 
one. A practical use is for sampling their sets, so getting bounds for a stable distance between 
size functions. Our extension will produce the same opportunity for persistent homology. 

2.2 1-dimensional Persistent Homology 

Given a topological space X and an integer i, we denote the i-th singular homology module of 
X over afield k by H,{X). 

Next we report two definitions of [5]. 

Definition 4. Let X he a topological space and f a real function on X. A homological critical 
value of f is a real number a for which there exists an integer i such that, for all sufficiently 
small e > 0, the map Hi{f ^^{—oo, a — e]) —>■ Hi{f ^^{—oo, a + e]) induced by inclusion is not an 
isomorphism. 

This is called an i-essential critical value in the paper [51 Def.2.6], dedicated to the size 
functor, a contemporary and not too different homological generalization of size functions. 

Definition 5. A function f : X ^ R is tame if it has a finite number of homological critical 
values and the homology modules Hi(f^^{^oo, a]) are finite- dimensional for alii Eli and a G M. 

The reader should be warned that there exist other, different meanings of "tame" in the 
current topological literature. Actually, "homologically tame" might be a better designation for 
such a type of function, but we adhere to this already current definition. 

We write i^" = Hi{f-^{-oo, u]), for all i e 1, and for u<v,we let : ^ be the 
map induced by inclusion of the lower level set of u in that of v, for a fixed integer i. Moreover, 
we indicate with F^^'" = Im/"'" the image of i^/' in F'" , that is called i-th persistent homology 
module. 
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3 Homological 1-dimensional reduction 



In this section we define the z-th persistent homology module related to a continuous n-dimensional 
real function (substantially as in [3J). Then we show that the sets of points of R^", where the 
modules change, can be obtained by computing the discontinuity points of persistent homology 
of a 1-dimensional function defined on particular half-planes which foliate the 2ri-space. 

The first issue arises when one tries to compute the maximum between the components of a 
n-dimensional real function. In fact: 

Remark 1. The maximum of two tame functions is not necessarily a tame function. 

(We recall that "tame" has the meaning defined in Section [^21) 
As an example, let /i, /2 : ^ M be two tame functions defined as 

— w — u^sin(i) u^O 

—V u = 

and consider the function 

/ = max(/i,/2). 

Then, as we can see in Figure [31 / is not tame, since iJo(/^^(— oo, 0]) is an infinite- 
dimensional module. 

Given this fault related to tame functions, the solution we propose is to introduce the fol- 
lowing concept. 

Definition 6. Let X be a topological space and f : X M" a continuous function on X. 

We shall say that f is max-tame if, for every admissible pair {l,b), the function g{P) = 

max \ i^lID. — h. I js tame. 
j=l,...,n I 'j J 

Choosing a measuring function on X as above, let us define the multidimensional persistent 
modules. 

Definition 7. Let f : X K" be a max-tame function. For each homology degree z G Z we 
put Ff = iJj(/-i(n (-oo,Uj])), for all u G R". For u ^ v we let : Ff Ff be the 
map induced by inclusion of the lower level set of u in that of v, for a fixed integer i, and call 
pu,v _ ju,v multidimensional persistent homology module. 

Note that the rank of F^'" is what is called px,i{u,v) in [31 Def. 12]. 



V — sa\{^) u ^ Q 



M = 



, /2(m,w) 
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II 







Figure 1: Lower level set of /i 
(grey area - one connected compo- 
nent) . 



Figure 2: Lower level set of /2 
(grey area - one connected compo- 
nent) . 



Figure 3: Lower level set of / (dark zone - infinitely many connected components). 
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Let g{P) = max j^^^^^^} for a fixed {l,b) £ Adm„, Gf = H,{g-^-oo, s]), for all s G M 

J — l,...,n L ^ J 

and « e Z. For s < t, we let g^'* : Gf — > G* be the map induced by inclusion of the lower level 
set of s in that of t, for a fixed integer i, and denote G^'* = Im g'^''^ the i-th persistent homology 
module. 

Now we can state and prove the theorem which, in analogy with the main result of [T], enables 
us to reduce the computation of multidimensional persistent homology to the 1-dimensional one. 
This is important, not so much for finding the homology modules themselves point by point, 
but much more for finding points of change of the modules. 

Theorem 2. Let (l,b) be an admissible pair and f — (/i,...,/„) : X M" a max-tame 
function. Then, for every (u, w) — {si + b,tl + b) £ ""(fS); following equality 

holds for all i Cz and s, t G M with s < t. 

Proof. By Lemma [1] we know that, for every u — (ui, . . . , Wn) S M", with uj — slj + bj,j — 
1, . . . , n, it holds that 

{P e X, f,{P) < u„j = 1, . . . , n} = {P e X,g{P) < s} 

hence 

{PeX,Pe ff\-^, J = 1, . . . , n} = {P e X, P e ^-^(-cx., s]}. 

It follows that 

n 

implying 

|^n/r'(-^,«j]j =H,{g-\-^,s\) 

for all i ^Ij. 

Analogously, for every v = {vi,...,Vn) G K", with Vj = tlj +bj,j = l,...,n, it holds 

/n _ \ _ 

that Hi Pi fj {—oo,Vj] = Hi{g oo,t]), for all i £ Z. So, since and g^' have the 

same domain and codomain and they are the maps induced by inclusion, we can conclude that 
— g*'*, and the claim follows. □ 
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4 Multidimensional matching distance 



According to [3, Def. 12], for a given measuring function f : X ^ M", for each homology degree 
i G Z the rank invariant ^ : — > N is defined as p'-^ ^{u, v) = rank{F^'^). 

Let {XJ'), {YJ") be two size pairs, where f':X-> R", /" : F -> R" are max-tame 
measuring functions, and "p'x T^'vi be the respective rank invariants. Let an admissible pair 
{f, b) be fixed, and let g' : X ^R, g" : y ^ R be defined by setting 

5 (P) = max <^ ^ g (P) = max <^ ^ 

] = l,...,n Ij J j = l,...,n J 

It is well-known for 1-dimensional measuring functions |181 113[ [5] that the relevant informa- 
tion on the rank invariants p'x i, Py i of g' and g" respectively is contained, for each degree i, 
in their multisets of cornerpoints, which are called "persistence diagrams" . These are sets of 
points of the extended plane with multiplicities, augmented by adding a countable infinity of 
points of the diagonal y ^ x: let them be called respectively C and C". Each cornerpoint is 
determined by its coordinates x < y < oo. The distance of two cornerpoints is 

S ((a, b), (c, d)) = min •|max{|a — c|, |6 — d\}, max 
It has been proved in [S] that the matching (or bottleneck) distance 

d{Px,i,PY,i) = min max (5(P,cr(P)) 
(T PeC 

where a varies among all bijections from C" to C", is stable. Mimicking [1] (and recalling that 
p'x iT Py i vary with (l,b)) we can use d to define distances between the rank invariants of the 
original multidimensional persistent homologies. 

Definition 8. Let {X, /'), (Y, /") be two size pairs and j, py ^ be the respective rank invari- 
ants. Then the i-th multidimensional matching distance between rank invariants is defined as 
the extended distance 

D{p'x,z,PY,i) ^ sup min Ij ■ d{p'x,„ Py^i) 

Note that D is by construction a global distance, i.e. not depending on {l,b), but since the 
coefficients Ij are < 1, there might be distances d, for particular admissible pairs, which take 
greater values. An easy corollary of our Theorem [2] is the following, which is the higher degree 
version of [Ii Cor. 1]. 



b — a d — 
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Corollary 1. For each z G Z the identity p'^i = p'y t holds if and only if d{p'j^ ^, py ^) = for 
every admissible pair {l,b). 

With the same argument of the analogous Proposition 4 of [IJ, it is easy to prove the foUowing 
inequahty between the muhidimensional matching distance and the 1-dimensional one obtained 
by considering the components of the measuring functions. That this inequahty can be strict, 
is shown in Section [5] 

Proposition 1. Let {X,f),{Y,h) be size two pairs with f = (/i....,/„), h = (/ii, . . . , /i„) 
max-tame measuring functions. For each i ^ TL and for each j = 1, . . . ,7i let Px iiPvi 
i-th rank invariants relative to the components fj, hj respectively; let then'p'xi, i be the rank 
invariants relative to /, h respectively. Then it holds that 

The matching distance is known to be stable with respect to perturbation of 1-dimensional 
measuring functions [5l Section 3.1] [7, Thm. 25]. In the multidimensional setting, the stability 
of d with respect to an admissible pair is stated in the following proposition, whose proof is 
again a copy of that of 1, Prop. 2]. Here 5', g" : X M. are defined in correspondence to (^, b) 
as at the beginning of this Section. 

Proposition 2. // {X,f'), {X,f") are size pairs, with max-tame functions f',f" : X E", 
and maxpgx ||/'(^') — f"iP)\\oo < then for every admissible pair (l, b) and for each i £1, it 
holds that 

mmj=i^...^„ Ij 

with I = (Zi, . . . , In) and where p'x j, p'x i the rank invariants at degree i of {X,g'), {X,g") 
respectively. 

By the definition of D, every 1-dimensional matching distance obtained in correspondence 
of an admissible pair yields a lower bound for the multidimensional matching distance D; a 
sufficiently fine sampling by admissible pairs produces approximations of arbitrary precision of 
it. 

Of course, maxD{'p'x i, Py J is still a meaningful distance related to the size pairs {X, /'), {Y, /"). 
We plan to study its relation with the n-dimensional natural pseudodistance of [14^ [1] . 
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5 Examples and Remarks 



We now describe a simple example, which shows that persistent homology, with respect to 
a mnlticiimensional measuring function, is actually stronger than the simple collection of the 
persistent homologies with respect to its l-climensional components. 

In M'^ consider the set O = [—1, 1] x [—1, 1] x [—1, 1] and the sphere § of equation w^ + y^+w;^ = 
1. Let also / = (/i, /2) : R'^ be a continuous function, defined as f{u, v, w) = {\u\, \v\). In 

this setting, consider the size pairs (C, ip) and (§, tp), where C = dil and tp and ijj are respectively 
the restrictions of / to C and S. 

In order to compare the persistent homology modules of C and § defined by /, we are 
interested in studying the half-planes' foliation of R^, where I = (cos0, sin^) with 6 e (0, |), 
and b = (a, —a) with a € R. Any such half-plane is parameterized as 

ui = s cos 6 + a 
U2 = ssinO — a 
vi = t cos 9 + a 
V2 = tsinO — a 

with s,t €R,s <t. 

In the following, we shall always assume < s <t. 

For example, by choosing d = j and a = 0, i.e. I = ^) and b = (0, 0), we obtain that 
g' = v^maxjiyji, (^2} = V^maxduj, \v\}, 
g" = \/2max{V'i, V'2} = V^maxjlul, 

Figure 4: Lower level sets g' < 1 and g" < I. 





11 



Let pg i p'li he the rank invariants of the respective persistent homologies for i G Z. So, 
writing G|(e) = Hi ((5')"^ (-00, s]), G|(S) = Hi ((5")"^ (-00, s]) and Gf defined as above, 
we obtain that ^ 

" 1 J 



0, 


s,i < 




< s < t < a/2 


k, 


otherwise 


0, 


s,t<0 




< s < i < 1 


fc, 


otherwise 



i=V2 





t=l 



t 



Gf(e)= 0, foralls,te 



t, 



(e) 



G^'*(S) 



A:^, 1 < s < t < 
0, otherwise 

k, \/2<s <t 
0, otherwise 

fc, V2 < .s < t 
0, otherwise 



t=V2 



2 



M 



S ^ i?(pe,2,ps,2) > ^d(Pe,2 Ps,2) = 



In other words, multidimensional persistent homology, with respect to and ip, is able to 
discriminate the cube and the sphere, while the 1-dimensional one, with respect to ipi , ^2 and 
tpi,tlJ2, cannot do that. In fact, for either manifold the lower level sets of the single components 
(i.e. 1-dimensional measuring functions) are homeomorphic for all values: they are topologically 
either circles, or annuli, or spheres. 

It should be noted that the map g" on S reaches the homological critical value 1 at points, 
at which it lacks of differentiability. 

In the example above, is not a Morse function (as would be desirable, if not necessary), 
because 5' : C — > M has infinitely many critical points when maxHul, \v\} = 1; moreover, the 
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cubic surface itself is not even C"'^. This problem can be solved by perturbing C so that it becomes 
smooth (e.g. a super-quadric |l6j). In this case, the differences between homology modules of 
the cube and of the super-quadric are only quantitative (i.e. the levels of homological critical 
values are different from one another). 

An even simpler example can be given on size pairs having the same support. Let X be 

u = cos 9 



^ + = 1 

the ellipse imbedded in M'^ as — or parameterized as < 



Let 



w = sin0 

(^1, (^2, V^i, : X ^ K. be defined as tpi = u, ipi = v, (p2 — ip2 — w and (f = {ipi,(p2), ip = 
(V'l; V'2)- Then the persistent homology modules of {X, ipi), {X, ijji), (X, ip2) — {X, 1P2) are iden- 
tical, while the persistent homology (in degree zero, so the size function) of (X, ip) differs from 
the one of (X, ijj). Indeed, while the lower level sets of ip are always either empty or connected. 



the lower level sets (/? < {u,w), with < u < 1, VT — < w < 1 consist of two connected 
components. 



6 Reduction of i-essential critical values 

The former example of the previous section suggests also some other considerations on the 
cooperation of measuring functions. We remind that the adjective "essential" is used here with 
the meaning introduced in 2^ and recalled in Section [2?2l after Definition HI so otherwise than, 
e.g., in g]. 

A first remark is that, although the persistent homology on single components of / cannot 
distinguish the two spaces, the persistent homology on /i restricted to lower level sets of /2 can, 
as can be shown as follows. Consider again the sphere §. The value l/\/2 (corresponding to the 
homological critical value 1 of g") is not critical for the maps /i, /2 on § itself, but it is indeed 
critical for /2 restricted to /j"^(— 00, l/\/2]. We believe that homological critical values of the 
ID reduction of multidimensional measuring functions are always clues of such phenomena. 

A further speculation on the use of cooperating measuring functions — from a completely 
different viewpoint than the one developed in the previous sections — is the following. A 
problem in 1-dimensional persistent homology, as well as for the size functor, is the computation 
of i-essential critical values for i > 0. A possibility is the use of several, independent measuring 
functions for lowering i, i.e. the degree at which the passage through the critical value causes 
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a homology change. Lowering i is important, since 0-essential critical values are easily detected 
by graph-theoretical techniques The following example shows that a suitable choice of a 
second, auxiliary measuring function may actually take 1-essential critical values to 0-essential 
ones. 

Let T be a torus of revolution around the x axis, with the innermost parallel circle of radius 
2, the outermost of radius 3. On 7 define (/i, /2) — {z, —z). Suppose we are interested in the 
persistent homology of the size pair (T, /i). Then (0, 0, 2) is a 1-essential critical point for /i, i.e. 
it is a point at which 1-degree homology changes. Of course, there are computational methods 
(e.g. by the Euler-Poincare characteristic) which enable us to detect it, but they will probably 
be tailored to the particular dimension of the manifold and to the particular homology degree. 

The same point is 0-essential for its restriction to /2^^(— oo, 1], so it can be recovered by the 
standard graph-theoretical techniques used in degree 0, i.e. for size functions. (The two functions 
need not be so strictly related: /2 could be replaced by Euclidean distance from (0, 0, 3) with 
the same effect). We conjecture that — at least whenever torsion is not involved — one can 
recursively take the i-essential values of a measuring function to (i — l)-essential ones, down to 
(easily computable) 0-essential critical values by means of other (auxiliary) measuring functions, 
as in this example. 

7 Conclusions and future work 

The need of extending persistent homology to the multidimensional case is a rather widespread 
belief, confirmed by simple examples. The present research shows the possibility of reducing the 
computation of persistent homology, with respect to multidimensional measuring functions, to 
the 1-dimensional case, following the line of thought of an analogous extension devised for size 
functions in [T. This reduction also yields a stable distance for the rank invariants of size pairs. 

In the next future, we plan to characterize the multidimensional max-tame measuring func- 
tions in a way that the reduction to ID case makes the specific features of persistent homology 
modules hold steady. It also would be our concern to give a rigorous definition of multidimen- 
sional homological critical values of a max-tame function and to relate them to the homological 
critical values of the maximum of its components. 

Eventually, in relation to our conjecture about i-essentiality (see SectionlH]), we plan to build 
an algorithm to recursively reduce j-essential critical points of a measuring function to 0-essential 
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ones. 
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